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J> ' Semileptonic and nonleptonic decays of the B c meson to charmonium and D 

mesons are studied in the framework of the relativistic quark model. The decay 
£f~^ . form factors are explicitly expressed through the overlap integrals of the meson wave 

functions in the whole accessible kinematical range. The relativistic meson wave 
functions are used for the calculation of the decay rates. The obtained results are 



compared with the predictions of other approaches. 
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43 I. INTRODUCTION 



The investigation of weak decays of mesons composed of a heavy quark and antiquark 
gives a very important insight in the heavy quark dynamics. The properties of the B c 
meson are of special interest, since it is the only heavy meson consisting of two heavy quarks 
with different flavor. This difference of quark flavors forbids annihilation into gluons. As 
a result, the excited B c meson states lying below the BD production threshold undergo 
pionic or radiative transitions to the pseudoscalar ground state which is considerably more 
stable than corresponding charmonium or bottomonium states and decays only weakly. The 
Collider Detector at Fermilab (CDF) Collaboration [1] reported the discovery of the B c 
ground state in pp collisions. More experimental data are expected to come in near future 
from the Tevatron and Large Hadron Collider (LHC). 

The characteristic feature of the B c meson is that both quarks forming it are heavy and 
thus their weak decays give comparable contributions to the total decay rate. Therefore it 
is necessary to consider both the b quark decays b — > c, u with the c quark being a spectator 
and c quark decays c —>■ s,d with b quark being a spectator. The former transitions lead 
to semileptonic decays to charmonium and D mesons while the latter lead to decays to 
B s and B mesons. The estimates of the B c decay rates indicate that the c quark decays 
give the dominant contribution (~ 70%) while the b quark decays and weak annihilation 
contribute about 20% and 10%, respectively (for a recent review see e.g. [2] and references 
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therein). However, from the experimental point of view the B c decays to charmonium are 
easier to identify. Indeed, CDF observed B c mesons [1] analysing their semileptonic decays 
B c -> J/iplu. 

The important difference between the B c semileptonic decays induced by b — > c, u and 
c — > s, d transitions lies in the substantial difference of their kinematical ranges. In the 
case of B c decays to charmonium and mesons the kinematical range (the square of 
momentum transfer to the lepton pair varies from to g^ ax ~ 10 GeV 2 for decays to J/ip 
and g^ ax 18 GeV 2 for decays to D mesons) is considerably broader than for decays to 
B^> and B^ mesons (g 2 iax ~ 0.8 GeV 2 for decays to B s and g 2 iax rs 1 GeV 2 for decays to 
B mesons). As a result in the B c meson rest frame the maximum recoil momentum of the 
final charmonium and D mesons is of the same order of magnitude as their masses, while 
the maximum recoil momentum of the B^ and B^ mesons is considerably smaller than 
the meson masses. This significant difference in kinematics makes it reasonable to consider 
B c decays induced by b and c quark decays separately. 

In this paper we consider weak B c decays to charmonium and D mesons in the framework 
of the relativistic quark model based on the quasipotential approach in quantum field theory. 
This model has been successfully applied for the calculations of mass spectra, radiative and 
weak decays of heavy quarkonia and heavy-light mesons [3-8]. In our recent paper [9] we 
applied this model for the investigation of properties of the B c meson and heavy quarkonia. 
The relativistic wave functions obtained there are used for the calculation of the transition 
matrix elements. The consistent theoretical description of B c decays to charmonium and D 
mesons requires a reliable determination of the q 2 dependence of the decay amplitudes in 
the whole kinematical range. In most previous calculations the corresponding decay form 
factors were determined only at one kinematical point either q 2 = or q 2 = g^ ax and then 
extrapolated to the allowed kinematical range using some phenomenological ansatz (mainly 
(di)pole or Gaussian). Our aim is to explicitly determine the q 2 dependence of form factors 
in the whole kinematical range in order to avoid extrapolations thus reducing uncertainties. 
The large values of recoil momentum require the consistent relativistic treatment of these 
decays. In particular, the relativistic transformation of the meson wave functions from the 
moving to the rest reference frame should be taken into account. On the other hand, the 
presence of only heavy quarks in B c and charmonium allows one to use expansions in the 
inverse powers of heavy quark masses l/m b c . 

The paper is organized as follows. In Sec. II we describe the underlying relativistic 
quark model. The method for calculating matrix elements of the weak current for b — > c, u 
transitions in B c meson decays is presented in Sec. III. Special attention is devoted to the 
dependence of the decay amplitudes on the momentum transfer. The B c decay form factors 
are calculated in the whole kinematical range in Sec. IV. The q 2 dependence of the form 
factors is explicitly determined. These form factors are used for the calculation of the B c 
semileptonic decay rates in Sec. V. Section VI contains our predictions for the energetic 
nonleptonic B c decays in the factorization approximation, and a comparison of our results 
with other theoretical calculations is presented. Our conclusions are given in Sec. VII. 
Finally, the Appendix contains complete expressions for the decay form factors. 

II. RELATIVISTIC QUARK MODEL 

In the quasipotential approach a meson is described by the wave function of the bound 
quark- ant iquark state, which satisfies the quasipotential equation [10] of the Schrodinger 
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type [11] 

(tr - a •««■>- I^v ipM ^ (d 

where the relativistic reduced mass is 

E,E 2 M 4 - (mj - m 2 ) 2 

and i?!, £" 2 are the center of mass energies on mass shell given by 

_ M 2 -ml + mj _ M 2 - m\ + m\ 

1_ 2M ' 2 ~ 2M ' ^ 

Here M = E\ + i? 2 is the meson mass, mi^ are the quark masses, and p is their relative 
momentum. In the center of mass system the relative momentum squared on mass shell 

2 = [M 2 - (mi + m 2 ) 2 ][M 2 - (m t - m 2 ) 2 ] 
1 J 4M 2 ' 1 J 

The kernel V(p, q; M) in Eq. (1) is the quasipotential operator of the quark- ant iquark in- 
teraction. It is constructed with the help of the off-mass-shell scattering amplitude, projected 
onto the positive energy states. Constructing the quasipotential of the quark-antiquark in- 
teraction, we have assumed that the effective interaction is the sum of the usual one-gluon 
exchange term with the mixture of long-range vector and scalar linear confining potentials, 
where the vector confining potential contains the Pauli interaction. The quasipotential is 
then defined by [3] 

V(p, q; M) = M 1 (p)M 2 (-p)V(p, q; M) Ul (q)u 2 (-q), (5) 

with 

V(p, q; M) = ^D^Hrt + (k)r$T 2;/1 + Ko„f (k), 
where a s is the QCD coupling constant, is the gluon propagator in the Coulomb gauge 

£°°(k) = -g, ^(k) = - ^f) , D" = D* = 0, (6) 



and k = p — q; 7 M and u(p) are the Dirac matrices and spinors 



e(p) + m 



\ HP) \ e (p)+m 



I \ A - (7) 



Here cr and x X are the Pauli matrices and spinors; e{p) = \/p 2 + m 2 . The effective long-range 
vector vertex is given by 

r,(k) = 7, + ^<v*"> (8) 

where k is the Pauli interaction constant characterizing the long-range anomalous chromo- 
magnetic moment of quarks. Vector and scalar confining potentials in the nonrelativistic 
limit reduce to 

V v (r) = (1 - e)Ar + B, 

V s (r) = eAr, (9) 
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reproducing 

V COQ{ (r) = V s (r) + Vv(r) = Ar + B, (10) 

where e is the mixing coefficient. 

The expression for the quasipotential of the heavy quarkonia, expanded in v 2 /c 2 without 
and with retardation corrections to the confining potential, can be found in Refs. [3] and 
[4, 9], respectively. The structure of the spin-dependent interaction is in agreement with 
the parameterization of Eichten and Feinberg [12]. The quasipotential for the heavy quark 
interaction with a light antiquark without employing the expansion in inverse powers of the 
light quark mass is given in Ref. [5]. All the parameters of our model like quark masses, 
parameters of the linear confining potential A and B, mixing coefficient e and anomalous 
chromomagnetic quark moment k are fixed from the analysis of heavy quarkonium masses 
[3] and radiative decays [6]. The quark masses m& = 4.88 GeV, m c = 1.55 GeV, m u ^ = 0.33 
GeV and the parameters of the linear potential A = 0.18 GeV 2 and B = —0.16 GeV have 
usual values of quark models. The value of the mixing coefficient of vector and scalar 
confining potentials e — — 1 has been determined from the consideration of the heavy quark 
expansion for the semileptonic B — > D decays [7] and charmonium radiative decays [6]. 
Finally, the universal Pauli interaction constant k — —1 has been fixed from the analysis 
of the fine splitting of heavy quarkonia 3 Pj- states [3] . Note that the long-range magnetic 
contribution to the potential in our model is proportional to (1 + k) and thus vanishes for 
the chosen value of k — — 1. It has been known for a long time that the correct reproduction 
of the spin-dependent part of the quark- antiquark interaction requires either assuming the 
scalar confinement or equivalently introducing the Pauli interaction with k — — 1 [3, 4, 13] 
in the vector confinement. 



III. MATRIX ELEMENTS OF THE ELECTRO WEAK CURRENT FOR b -» c, u 

TRANSITIONS 

In order to calculate the exclusive semileptonic decay rate of the B c meson, it is necessary 
to determine the corresponding matrix element of the weak current between meson states. 
In the quasipotential approach, the matrix element of the weak current = ^7^(1 — 75)6, 
associated with b — > q (q = c or u) transition, between a B c meson with mass M Bc and 
momentum p Bc and a final meson F (F = ip,r) c or D^) with mass Mp and momentum pp 
takes the form [14] 

(F{p F )\J?\B e {p Be )) = J ^^^ P ,(p)r M (p,q)* BcPBc (q), (11) 

where T /1 (p,q) is the two-particle vertex function and $« PM are the meson (M = B C ,F) 
wave functions projected onto the positive energy states of quarks and boosted to the moving 
reference frame with momentum pm- 

The contributions to T come from Figs. 1 and 2. The contribution is the consequence 
of the projection onto the positive-energy states. Note that the form of the relativistic 
corrections resulting from the vertex function is explicitly dependent on the Lorentz 
structure of the quark- antiquark interaction. In the leading order of the v 2 /c 2 expansion 
for B c and tp and in the heavy quark limit m c — > 00 for D only T* 1 ) contributes, while r^ 2 ) 
contributes already at the subleading order. The vertex functions look like 

r^(P, q) = u q {p q )^(l - 7 5 Mg b )(27T) 3 5(p c - q c ), (12) 
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FIG. 1: Lowest order vertex function I^ 1 ) contributing to the current matrix element (11). 
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FIG. 2: Vertex function taking the quark interaction into account. Dashed lines correspond to 
the effective potential V in (5). Bold lines denote the negative-energy part of the quark propagator. 



and 



rf (p,q) = u q (p q )u c ( Pc ){^(l - 7l ) 



e b (k) + e b (p q ) 



7i°V( Pc - q c ) 



+V( Pc - q c ) 



e q (k') +e q (q b ) 



TiTnX 1 - li)} u b(qb)u c (q c ), (13) 



where the superscripts "(1)" and "(2)" correspond to Figs. 1 and 2, k = p q — A; k' 
q b + A; A = Pf — p Bc ; 



A(")(p) 



e(p) - (m7° + 7 (7p)) 

Hp) 



Here [14] 



Pq,< 



Qb,, 



F i=l 

zUq)^±i:^(PB C )q\ 



B c i=i 



and are three four-vectors given by 



M 



M(E + M) 



j, E = \Jp 2 + ~AP. 
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It is important to note that the wave functions entering the weak current matrix element 
(11) are not in the rest frame in general. For example, in the B c meson rest frame (p# c = 0), 
the final meson is moving with the recoil momentum A. The wave function of the moving 
meson f f A is connected with the wave function in the rest frame ^fo = by the 
transformation [14] 

* fa(p) = (it>rjz^wj^o( P ), (M) 

where R w is the Wigner rotation, La is the Lorentz boost from the meson rest frame to a 
moving one, and the rotation matrix D l l 2 (R) in spinor representation is given by 



is the usual Lorentz transformation matrix of the four-spinor. 

The general structure of the current matrix element (11) is rather complicated, because it 
is necessary to integrate both with respect to d 3 p and d 3 q. The 5-function in the expression 
(12) for the vertex function I^ 1 ) permits to perform one of these integrations. As a result 
the contribution of I^ 1 ) to the current matrix element has the usual structure of an overlap 
integral of meson wave functions and can be calculated exactly (without employing any 
expansion) in the whole kinematical range, if the wave functions of the initial and final 
meson are known. The situation with the contribution is different. Here, instead of a 
5-function, we have a complicated structure, containing the potential of the gg-interaction 
in meson. Thus in the general case we cannot get rid of one of the integrations in the 
contribution of r^ 2 ^ to the matrix element (11). Therefore, it is necessary to use some 
additional considerations in order to simplify calculations. The main idea is to expand 
the vertex function T^, given by (13), in such a way that it will be possible to use the 
quasipotential equation (1) in order to perform one of the integrations in the current matrix 
element (11). 

B c — > if), r) c transitions. The natural expansion parameters for B c decays to charmo- 
nium are the active heavy b and c quark masses as well as the spectator c quark mass. We 
carry such an expansion up to the second order in the ratios of the relative quark momentum 
p and binding energy to the heavy quark masses m^ c . It is important to take into account the 
fact that in the case of weak B c decays caused by b — > c, u quark transition the kinematically 
allowed range is large (|A max | = [M\ c — M f )/(2Mb c ) ~ 2.4 GeV for decays to charmonium 
and ~ 2.8 GeV for decays to D mesons). This means that the recoil momentum A of a 
final meson is large in comparison to the relative momentum p of quarks inside a meson 
(~ 0.5 GeV), being of the same order as the heavy quark mass almost in the whole kinemati- 
cal range. Thus we do not use expansions in powers of | A|/m(, )C or | A|/Mp, but approximate 

in the expression (13) for the heavy quark energies e^ c (p + A) = \Jml c + (p + A) 2 by 



ef, )C (A) = y m 2 c + A 2 , which become independent of the quark relative momentum p. Mak- 
ing these replacements and expansions we see that it is possible to integrate the current 
matrix element (11) either with respect to d 3 p or d 3 q using the quasipotential equation (1). 
Performing integrations and taking the sum of the contributions of and we get 
the expression for the current matrix element, which contains ordinary overlap integrals of 




(15) 



where 
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FIG. 3: The reduced radial wave functions for the 5-states of the B c meson. 

meson wave functions and is valid in the whole kinematical range. Thus this matrix element 
can be easily calculated using numerical wave functions found in our meson mass spectrum 
analysis [4, 9]. The reduced radial wave functions u(r) = rR(r) of the B c meson are shown 
in Fig. 3. 

B c — ► transitions. In this case the heavy b quark undergoes the weak transition 
to the light u quark. The constituent u quark mass is of the same order of magnitude as 
the relative momentum and binding energy, thus we cannot apply the expansion in inverse 
powers of its mass. Nevertheless, taking into account the fact that the recoil momentum 
of the final meson in this decay is large almost in the whole kinematical range (as it was 
discussed above), we can neglect the relative momentum p of quarks inside a meson with 
respect to the large recoil momentum A. Thus in the region of large recoil (|A| ^> |p|) 
we can use the same expressions of the r^ 2 ^ contribution to the current matrix element 
both for the B c — > and B c — > ^,r] c transitions. Moreover, the smallness of the r^ 2 ^ 
contribution, which is proportional to the small binding energy, and its weak dependence 
on momentum transfer allows one to extrapolate these formulae to the whole kinematical 
range. As numerical estimates show (see below), such extrapolation introduces only small 
uncertainties. 



IV. B r DECAY FORM FACTORS 



The matrix elements of the weak current J w for B c decays to pseudoscalar mesons (P 
rjc, D) can be parametrized by two invariant form factors: 



(P(p F )\<rfb\B c (p BB )) = f + (q 2 ) 



' „ M 2 B - Ml " 

Pk+PF- Bc q2 P <t 



M 2 — M 2 

+ /o(g 2 ) MB < q2 Mj V , (16) 



where q = Pb c ~ Pf] Mb c is the B c meson mass and Mp is the pseudoscalar meson mass. 
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The corresponding matrix elements for B c decays to vector mesons (V = J/ip,D*) are 
parametrized by four form factors 

(V(p F )\qYb\B(p Bc )) = M^fl v ^ pa <PB. P PF«, (17) 

(V( PF )\qri 5 b\B(p Bc )} = 2M v A (q 2 ) € —^ # + (M Bc + M v )A l (q 2 ) (e*» - 



q 

-A 2 {q 2 ) 




Mr,. + M v 



u Ml - My , 

p% c +pf- c 2 v q" 

q z 



(18) 



where My and e M are the mass and polarization vector of the final vector meson. The 
following relations hold for the form factors at the maximum recoil point of the final meson 
{q 2 = 0) 

/+(0) = /o(0), 

In the limit of vanishing lepton mass, the form factors /o and A do not contribute to 
the semileptonic decay rates. However, they contribute to nonleptonic decay rates in the 
factorization approximation. 

It is convenient to consider B c semileptonic and nonleptonic decays in the B c meson rest 
frame. Then it is important to take into account the boost of the final meson wave function 
from the rest reference frame to the moving one with the recoil momentum A, given by 
Eq. (14). Now we can apply the method for calculating decay matrix elements described 
in the previous section. As it is argued above, the leading contributions arising from the 
vertex function can be exactly expressed through the overlap integrals of the meson 
wave functions in the whole kinematical range. For the subleading contribution F^ 2 \ the 
expansion in powers of the ratio of the relative quark momentum p to heavy quark masses 
ffli, >c should be performed taking into account that the recoil momentum of the final meson 
A can be large. Such expansion is well justified for B c decays to charmonium in the whole 
kinematical range. For B c decays to D mesons, where one of the final quarks is light, a 
similar expansion is well justified only in the kinematical region of large recoil momentum. 
However, the numerical smallness of this subleading contribution due to its proportionality 
to the small meson binding energy permits its extrapolation to the whole kinematical range. 
As a result, we get the following expressions for the B c decay form factors: 

(a) B c — > P transitions (P = i] c , D) 

U(q 2 ) = ff{q 2 )+ef S + {2 \q 2 ) + (1 - e )/£<V), (19) 

fo(q 2 ) = / (1) (g 2 ) + ef S(2 \q 2 ) + (1 - £)/ y( V), (20) 

(b) B c -> V transition (V = ip, D*) 

V (q 2 ) = V (1 \q 2 ) + eV S(2 \q 2 ) + (1 - e)V V(2 \q 2 ), (21) 

Mq 2 ) = A^(q 2 ) + eA^ 2 \q 2 ) + (1 - e)A V ^{q% (22) 
A,{q 2 ) = A?(q 2 )+eA^ 2 \q 2 ) + (1 - e)A v ^ 2 \q\ (23) 
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TABLE I: Form factors of weak B c decays (6 — > c,u transitions). 



Transition 




f < «2\ 

M<? J 




A ( „2\ 




A o(q ) 


/ ) . __ T /„/, 














9 9max 


l.U < 


n no 
U.yz 


l.o4 


U.oo 


1.66 


i.uo 


g 2 = 


0.47 


0.47 


0.49 


0.50 


0.73 


0.40 
















Q (/max 


0.08 


0.05 


-0.16 


0.03 


0.10 


0.08 


g 2 = 


0.27 


0.27 


0.24 


0.18 


0.14 


0.23 


— ► D, D* 
















1.20 


0.64 


2.60 


0.62 


1.78 


0.97 


q 2 =0 


0.14 


0.14 


0.18 


0.17 


0.19 


0.14 




A (q 2 ) 


= 4V) + ^ 5( V) + 


(l-e)< ( V), 




(24) 



where /+o, f+o , ^012? A) '12 ■> an d V s ' v ^ are given in Appendix. The superscripts 
"(1)" and "(2)" correspond to Figs. 1 and 2, S and V to the scalar and vector potentials of 
gg-interaction. The mixing parameter of scalar and vector confining potentials e is fixed to 
be —1 in our model. 

It is easy to check that in the heavy quark limit the decay matrix elements (16)-(18) 
with form factors (19)-(24) satisfy the heavy quark spin symmetry relations [15] obtained 
near the zero recoil point (A — > 0). 

For numerical calculations we use the quasipotential wave functions of the B c meson, 
charmonium and D mesons obtained in the mass spectra calculations [4, 5]. Our model 
predicts the B c meson mass Mb c = 6.270 GeV [9], while for J/ip, f] c , ip', rf c , D and D* 
meson masses we use experimental data [16]. The calculated values of form factors at zero 
(q 2 = g^ax) an d maximum (q 2 = 0) recoil of the final meson are listed in Table I. In Fig. 4 
we plot leading and subleading V s ^ 2 \ V v ^ contributions to the form factor V for 
B c — > D* transition, as an example. We see that the leading contribution is dominant 
in the whole kinematical range, as it was expected. The subleading contributions V s ^ 2 \ 
V v ^ are small and depend weakly on q 2 . The behavior of corresponding contributions to 
other form factors is similar. This supports our conjecture that the formulae (A.1)-(A.18) 
can be applied for the calculation of the form factors of B c — > transitions in the whole 
kinematical range. 

In Figs. 5-10 we plot the calculated q 2 dependence of the weak form factors of Cabibbo- 
Kobayashi-Maskawa (CKM) favored B c — > r] c , B c — > J/ip, B c — > rj' c , B c — > ip', as well 
as CKM suppressed B c — > D, B c — > D* transitions in the whole kinematical range. The 
different behavior (growing or falling with q 2 ) of the form factors displayed in Figs. 5-8 is 
evoked by the properties of the final meson wave functions, since the 2S wave function of 
the radially excited 77', ip' mesons has a zero. 

In the following sections we use the obtained form factors for the calculation of the B c 
semileptonic and nonleptonic decay rates. 
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FIG. 4: Leading V^' and subleading V s ^ 2 \ V v ^ contributions to the form factor V for the 
B c —* D* transition. 
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FIG. 5: Form factors of B c — > ij c eu decay. 

V. SEMILEPTONIC DECAYS 

The differential semileptonic decay rates can be expressed in terms of the form factors as 
follows. 
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FIG. 6: Form factors of B c — > J /ipev decay. 
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FIG. 7: Form factors of B c — > r^ez/ decay. 

(a) P c — > Pez/ decays (P = r/c, P) 

_ p e .) = (25 ) 

(b) P c -> Vez/ decays (V = ip, D*) 

- Vev) = (l» + (« 2 )| 2 + \H-tf)\ 2 + Iffol, 2 )! 2 ) , (26) 
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FIG. 8: Form factors of B c 
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FIG. 9: Form factors of B c — * Dev decay, 
where Gp is the Fermi constant, V q b is the CKM matrix element (q = c,u), 

A = IAI 



(M 2 Bc + My-q 



2\2 



The helicity amplitudes are given by 

w 2Mb - a 



4M| 



F(?2)T (M flc + M y )^^ 



2Mb A 



-M<f 



(27) 
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H (q 2 ) 1 



AM 2 A 2 

(M Bc + M v )(M 2 Bc -M 2 V - q 2 )A x {q 2 ) ~ Mg + My M^ 



(28) 



The decay rates to the longitudinally and transversely polarized vector mesons are defined 
by 

dr L _ G 2 F ^\V qh \ 2 q 2 



dq 2 96vr 3 M\ 

dr T dr+ dr_ G\^\v qh \ 2 q 2 



H (q 2 )\\ (29) 



d , 2 dg2 + ^-^^4^ 2)|2+l ^ (g2)|2 )- (30) 

In Figs. 11-16 we plot the differential semileptonic decay rates dT/dq 2 for semileptonic 
decays B c — > rj c (J /i/j)ev, B c — > j]' c {ijj')ev and B c — > D(D*)ei' calculated in our model using 
Eqs. (25), (26) both with and without account of l/m^ c corrections to the decay form 
factors (A.1)-(A.18). 1 From these plots we see that relativistic effects related to heavy 
quarks increase semileptonic B c decay rates to the pseudoscalar T] c , rj' c and D mesons, while 
semileptonic decay rates to vector J/ip, ^' and D* mesons are decreased by them. The 
decay rates for B c — > v( c ev and B c — ► Dev receive the largest l/m^^ corrections. This is not 
surprising since in the former decay the radially excited r)' c wave function has a zero, which 
considerably decreases the nonrelativistic contribution and thus increases the relative size 
of relativistic effects. In the latter decay, the role of relativistic effects is enhanced due to 
the relativistic light quark in the D meson. 

We calculate the total rates of the semileptonic B c decays to the ground and radially 
excited states of charmonium and D mesons integrating the corresponding differential decay 
rates over q 2 . For calculations we use the following values of the CKM matrix elements: 
\Vcb\ = 0.041, \ V u b\ = 0.0036. The results are given in Table II in comparison with predictions 



Relativistic wave functions were used for both calculations. 
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FIG. 11: Differential decay rates (l/|Kt>| 2 )dr/dg 2 of B c -» rj c ev decay (in 1(T 12 GeV" 1 ). The 
lower curve is calculated without account of l/m(, )C corrections. 
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FIG. 12: Differential decay rates (l/|y cil | 2 )dr/d<? 2 of B c J/^eu decay (in 1(T 12 GeV" 1 ). The 
upper curve is calculated without account of l/mfe jC corrections. 

of other approaches based on quark models [17, 19, 20, 22, 23, 25], QCD sum rules [18] and 
on the application of heavy quark symmetry relations [21, 24] to the quark model. Our 
predictions for the CKM favored semileptonic B c decays to charmonium ground states are 
almost 2 times smaller than those of QCD sum rules [18] and quark models [17, 19, 20], but 
agree with quark model results [22-25] . Note that the ratios of the B c — > J / ipeu to B c — ► r\ c ev 
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FIG. 13: Differential decay rates (l/\V cb \ 2 )dT/dq 2 of B c -» rf c ev decay (in 1(T 12 GeV- 1 ). The 
lower curve is calculated without account of l/m(, )C corrections. 
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FIG. 14: Differential decay rates {l/\V cb \ 2 )dT /dq 2 of B c tf/eu decay (in 1(T 12 GeV" 1 ). The 
upper curve is calculated without account of l/m^ c corrections. 

decay rates have close values in all approaches except [21]. In the case of semileptonic decays 
to radially excited charmonium states our prediction for the decay to the pseudoscalar rj c 
state is consistent with others, while the one for the decay to ip' is considerably smaller (with 
the exception of Ref. [25]). For the CKM suppressed semileptonic decays of B c to D mesons 
our results are in agreement with those of Ref. [21]. 
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FIG. 15: Differential decay rate (l/\V ub \ 2 )dT/dq 2 of B c -> Dev decay (in 1(T 12 GeV" 1 ). The lower 
curve is calculated without account of l/m& iC corrections. 
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FIG. 16: Differential decay rates (l/|K6| 2 )dr/dg 2 of B c -> D*ei/ decay (in 1(T 12 GeV" 1 ). The 
upper curve is calculated without account of l/mb jC corrections. 



In Table III we present for completeness our predictions for the rates of the semileptonic 
B c decays to vector {ip and D*) mesons with longitudinal (L) or transverse (T) polarization 
and to the states with helicities A = ±1, as well as their ratios. 



17 



TABLE II: Semileptonic decay rates T (in 10 GeV) of B c decays to charmonium and D mesons. 

Decay our [17] [18] [19] [20] [21] [22] [23] [24] [25] 

B c -► r] c eu 5.9 14 11 11.1 14.2 2.1(6.9) 8.6 6.8 4.3 8.31 

B c ^r]' c ev 0.46 0.60 0.73 0.3 0.605 

B c -► J/ipev 17.7 33 28 30.2 34.4 21.6(48.3) 17.5 19.4 16.8 20.3 

B c -► ^ei/ 0.44 1.94 1.45 1.7 0.186 

B c ^Dev 0.019 0.26 0.059 0.049 0.094 0.005(0.03) 0.001 0.0853 

B c ^D*eu 0.11 0.49 0.27 0.192 0.269 0.12(0.5) 0.06 0.204 



TABLE III: Semileptonic decay rates 1^ ^,+,- (in 10 15 GeV) and their ratios for B c decays to 
vector tb and D* mesons. 



Decay 




r T 




r+ 


r_ 


r + /r_ 


B c — ► J/ibev 


7.8 


9.9 


0.78 


2.9 


7.0 


0.40 


B c -> 


0.29 


0.15 


1.85 


0.05 


0.10 


0.47 




0.04 


0.07 


0.53 


0.015 


0.055 


0.24 



VI. NONLEPTONIC DECAYS 



In the standard model nonleptonic B c decays are described by the effective Hamiltonian, 
obtained by integrating out the heavy VF-boson and top quark. For the case of b — > c, u 
transitions, one gets 



Cl (//)Of + c 2 {fj,)0, 



cb 



Cl (n)Of + c 2 (fi)O l 2 



+ .... 



(31) 



The Wilson coefficients Ci^/i) are evaluated perturbatively at the W scale and then are 
evolved down to the renormalization scale /x ~ by the renormalizat ion-group equations. 
The ellipsis denote the penguin operators, the Wilson coefficients of which are numerically 
much smaller than C\ 2 - The local four-quark operators 0\ and 2 are given by 



Of = [(du) V - A + (sc) V - A ](qb)v-A, 

Of = (qu) V - A (db)v„ A + (qc) v „ A (sb) V -A, q = (u,c), 



(32) 



where the rotated antiquark fields are 

d = V ud d + V us s, s = V cd d + V cs s, (33) 
and for the hadronic current the following notation is used 

{qq')v-A = ?7m(! ~ 7s)?' = J™ ' ■ 

The factorization approach, which is extensively used for the calculation of two-body 
nonleptonic decays, such as B c — > FM, assumes that the nonleptonic decay amplitude 
reduces to the product of a form factor and a decay constant [26]. This assumption in 
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B+ F° = tP, Vc 



FIG. 17: Quark diagram for the nonleptonic B+ — ► F°M + decay. 

general cannot be exact. However, it is expected that factorization can hold for the energetic 
decays, where the final F meson is heavy and the M meson is light [27]. A justification of 
this assumption is usually based on the issue of color transparency [28] . In these decays the 
final hadrons are produced in the form of point-like color-singlet objects with a large relative 
momentum. And thus the hadronization of the decay products occurs after they are too 
far away for strongly interacting with each other. That provides the possibility to avoid the 
final state interaction. A more general treatment of factorization is given in Refs. [29, 30]. 

In this paper we limit our analysis of the B+ nonleptonic decays to the case when the final 
meson F° is charmonium 2 and the light M + meson is ir + , p + or K^ + . The corresponding 
diagram is shown in Fig. 17, where q\ — d, s and q 2 = u. Then the decay amplitude can be 
approximated by the product of one-particle matrix elements 

(F°M + \H cS \B+) = ^V cb V qiq MF\(Hv-A\B c )(M\(q iq2 ) v _ A \0), (34) 



where 



and N r is the number of colors. 



ai = ci(/x) + — c 2 (//) (35) 



The matrix element of the current Jj/ between vacuum and final pseudoscalar (P) or 
vector {V) meson is parametrized by the decay constants fpy 

(^I?i7^75?2|0) = ifpft, (V\qa,q 2 \0} = e,M v f v . (36) 

We use the following values of the decay constants: f w = 0.131 GeV, f p = 0.208 GeV, 
f K = 0.160 GeV and f K . = 0.214 GeV. The CKM matrix elements are \V ud \ = 0.975, 
\V US \ = 0.222. 

The matrix elements of the weak current between the B c meson and the final charmonium 
entering the factorized nonleptonic decay amplitude (34) are parametrized by the set of decay 
form factors defined in Eqs. (16) and (17). Using the form factor values calculated in Sec. IV, 
we get predictions for the nonleptonic B+ — > F°M + decay rates and give them in Table IV 
in comparison with other calculations [18-22, 25]. We see that for most decays our model 
predicts slightly lower decay rates than other approaches. 



2 We do not consider nonleptonic B c decays where the final meson F is a D meson, since such decays 
are strongly CKM suppressed and thus receive important contributions from the weak annihilation and 
penguins. 
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TABLE IV: Nonleptonic B c decay rates T (in 1CT 15 GeV). 



Decay 


our 


[18] [19] 


[20] 


r^-»r\l 

[22] 


[21] 


[25] 


Bt 


-> 7] C TT + 


0.93a? 


1.8a? 1.59a? 


2.07a? 


1.47a? 


0.28a? 


1.49a? 


Bt 




2.3a? 


4.5a? 3.74a? 


r- A o 9 

5.48a? 


3.35a? 


0.75a? 


3.93a? 


B+ -► J/-07T+ 


0.67a? 


1.43a? 1.22a? 


1.97a? 


0.82a? 


1.48a? 


1.01a? 


B + -+ J/^ p + 


1.8a? 


4.37a? 3.48a? 


5.95a? 


2.32a? 


4.14a? 


3.25a? 




— ► r/c-K 


0.073a? 


0.15a? 0.119a? 


0.161a? 


0.15a? 


0.023a? 


0.115a? 


B+- 




0.12a? 


0.22a? 0.200a? 


0.286a? 


0.24a? 


0.041a? 


0.198a? 


B+- 


J/tpK + 


0.052a? 


0.12a? 0.090 a? 


0.152a? 


0.079a? 


0.076a? 


0.0764a? 


Bt~ 


> J/^K* + 


0.11a? 


0.25a? 0.197a? 


0.324a? 


0.18a? 


0.23a? 


0.174a? 


B+ 


-► f^7T + 


0.19a? 




0.268a? 




0.074a? 


0.248a? 


Bt 




0.40a? 




0.622a? 




r\ t n 9 

0.16a? 


0.587a? 


Bt 




0.12a? 




0.252a? 




0.22a? 


0.0708a? 


Bt 




0.20a? 




0.710 a? 




0.54a? 


0.183a? 


Bt 




0.014a? 




0.020a? 




0.0055a? 


0.0184a? 


Bt' 




0.021a? 




0.031a? 




0.008a? 


0.0283a? 


Bt 


->■ 


0.009a? 




0.018a? 




0.01a? 


0.00499a? 


Bt- 


-»• 


0.011a? 




0.038 a? 




0.03a? 


0.00909a? 



VII. CONCLUSIONS 

In this paper we calculated weak semileptonic and nonleptonic B c decays to the char- 
monium and D meson final states. The corresponding decay form factors were calculated 
in the framework of the relativistic quark model using v/c expansion for the B c meson and 
charmonium and heavy quark expansion for D mesons. These transitions proceed in a large 
kinematically allowed region. As a result, the recoil momentum A of the final meson in the 
B c rest frame is mostly large compared to the binding energy and the relative momentum 
of quarks forming a meson. Our approach permits to determine explicitly the dependence 
of form factors of the CKM favored B c transition to charmonium in the whole kinematical 
region through the overlap integrals of the meson wave functions. This is a real achievement 
making our results more reliable, since most other approaches determine form factors only at 
the single point of zero (q 2 = g^ ax ) or maximum (q 2 = 0) recoil of the final meson and then 
the extrapolation is used. We calculated form factors of weak B c transitions to charmonium 
up to second order corrections in the ratio of the relative quark momentum to the heavy 
(b and c) quark mass. In the case of the CKM suppressed B c transitions to D mesons the 
situation is more complicated, since the active quark undergoes heavy-to-light transition. 
The leading contribution as in the previous case can be determined exactly, while the sub- 
leading contribution, which is suppressed by the small binding energy, can be determined 
reliably in most part of the kinematical range except a small region near the point of zero 
recoil (g 2 = q 2 iax ) of the final D meson. As the numerical analysis shows, the extrapolation 
of form factors obtained in such a way to the region of small recoil introduces only minor 
uncertainties. 

We calculated semileptonic and nonleptonic (in factorization approximation) B c decay 
rates. Our predictions for the branching fractions are summarized in Table V, where we use 
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TABLE V: Branching fractions (in %) of exclusive B c decays calculated for the fixed values of the 
B c lifetime t Bc = 0.46 ps and a\ = 1.14. 



Decay 


Br 


Decay 


Br 


Decay 


Br 


B c ->■ r] c ev 


0.42 


B+ -► r/ c 7r+ 


0.085 




— ► ?/7r + 


0.017 


B c — >■ r]' c ev 


0.032 


B+ r/ c/ 9+ 


0.21 




->t/cP + 


0.036 


B c — ► J/ipev 


1.23 


5+ -► J/Y>7r+ 


0.061 


5 C + 




0.011 


B c — ► ip'ev 


0.031 


s+ - J/i>p + 


0.16 






0.018 


B c -► Dei/ 


0.013 




0.007 






0.001 


5c -► D*ezv 


0.037 


B+ _> r] c K* + 


0.011 






0.002 






B+ - J/ipK + 


0.005 




-► 


0.001 






B+ -► J/i>K* + 


0.010 


5+ -► v'if * + 


0.001 



the central experimental value of the B c meson lifetime [16]. From this table we see that 
the considered semileptonic decays to charmonium and D mesons give in total 1.72% of the 
B c decay rate, while the energetic nonleptonic decays give additional 0.63%. It is expected 
that the dominant contribution to the B c total rate comes from the charmed quark decays. 
These decays will be considered in a forthcoming publication. 
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APPENDIX: FORM FACTORS OF WEAK B c DECAYS 



(a) B c — > P transition (P = r) c , D) 



\M Bc 
x|l + 

X 



d3p <M P + 



2m r 



(2tt) 3 
M Bc - E P 
e q (p + A) + m q 
1 

+ 



e q (p + A) +m q 



+ 



Ep + Mp 2e q (p + A) \ 2e b (p) 

<pA) ( A 2 



e b (p) +m b 



e q (p + A)+m q e b (p) 
E P - M P 



—)) 



[e q (p + A) + m q ] [e b (p) + m h ] 



+ (M Bc - E P ) 



+ 



[e q {p + A) + m q ][e b (p) + m b ] 
1 1 N 



3 * \[tq(p + A) + m q ][e b (p) + m b ] \e q (p + A) + m q e c (p)+m c 
M Bc -E P ( 1 1 \ 



x 



Ep + Mp \e q (p + A)+m q e b (p) + m b J 
{e c (p)+m c ~ e 9 (p + A)+mJ)}* Bc(p) ' 



(A.l) 
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/ + 5(2 V) 



M, 



E P r d 3 p = / 2m c \ e g (A)+mJ 1 / e g (A) - 

(2tt) 3 P r^E P + Mp ;\ 2e,(A) \ e,(A) \e q (A) + 



M Bc -E P \ 
e q (A) +mq) 



Mi 



p + 



2m c 



(PA) 
A 2 



1 
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1 /e,(A)-m, 
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" ; Ve6(A)+m 5 e,(A) + m,j ' r " Dc " V e &( A ) + 

2m r 
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(b) B c -> V transition (V = V, £>*) 
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